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Abstract 



Positive powers of the mass parameter in a physical quantity calculated with 
the help of heavy quark effective theory originate from a Wilson coefficient in 
the matching of QCD and HQET Green function. We show that this mass 
parameter enters the calculation as a well-defined running current mass. We 
further argue that the recently found ill-definition of the pole mass, which is 
the natural expansion parameter of HQET, does not affect a phenomenological 
analysis which uses truncated perturbative series. We reanalyse inclusive 
semileptonic decays of heavy mesons and obtain the c quark mass m^ IS (m c ) = 
(1.35 ±0.20) GeV where the error is almost entirely due to scale-uncertainties. 
We also obtain mf(mj) = (4.6±0.3) GeV and \V cb \{T B /1.4Qy>s) 1/2 = 0.036± 
0.005 where the errors come from the uncertainty in the kinetic energy of the 
heavy quark inside the meson, in the experimental branching ratios, in QCD 
input parameters, and scale-uncertainties. 
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Typeset using REVIp^X 



I. INTRODUCTION 



During recent years the study of HQET, the effective theory of QCD expanded in inverse 
powers of the heavy quark mass, has considerably enlarged our understanding of low-energy 
QCD (cf. for reviews). Being successfully applied to a number of exclusive decay processes 
of heavy-light systems containing one heavy quark and one or two light quarks, it likewise 
supplies us with a number of relations between static properties of such systems like the 
particle spectrum and leptonic decay constants. Nevertheless, the question of the correct 
definition of the heavy quark mass tuq, the expansion parameter of HQET, has played only a 
minor role in the tremendous number of publications dealing with that theory. When Green 
functions in full QCD are matched to those in HQET, one has to decide whether to identify 
itlq with the pole mass or with the renormalization scheme-dependent current mass of the 
QCD Lagrangian. Admittedly the question of which mass to use seems a rather academical 
one from the viewpoint of phenomenology as long as the mass enters explicitly only in inverse 
powers and implicitly in the strong coupling constant in the matching coefficients. Yet the 
field of applications of HQET has expanded again and now likewise encloses semileptonic 
inclusive decays of heavy mesons @-fl- Two important statements about inclusive decays 
could be obtained: first, in leading order they are essentially free quark decays, and second, 
the leading corrections to the free quark decay are of order I/tuq. These results stimulated 
new determinations of the quark masses m c and and the CKM matrix element \V c b\ 
from the experimental measurements of the semileptonic branching ratios B(D — > Xev) 
and B(B — > X c eu), cf. Refs. 0-0]. Since in contrast to previous applications of HQET, 
the inclusive decay rate depends on the heavy quark mass in its fifth (!) power, a careful 
analysis of the correct definition of the mass parameter is mandatory. In fact any positive 
power of the heavy mass originates from the matching of a QCD Green function to a HQET 
Green function, so it is naturally a running current mass which enters in this step. On the 
other hand the expansion parameter of the HQET Green function equals the pole mass [|5],[| , 
whose precise definition has recently been investigated by Braun and Beneke and by Bigi et 
al. ||1U|| . It was found that due to renormalon effects any attempt to define the pole mass 
beyond a finite order in perturbation theory is plagued by an intrinsic uncertainty of order 
Aqcd- 

The present paper is organized as follows: In Sec. II we analyse the role of the mass 
parameter in the Wilson coefficient obtained from the matching of QCD and HQET Green 
functions. As a result we recover in the leading order of the 1 / mQ-expansion the long-known 
renormalization scheme- and scale-dependence of the parton model. We further argue that 
the ambiguity in the definition of the pole mass, which is the expansion parameter of HQET, 
does not affect a phenomenological analysis which uses the expansions in 1/ uiq and a s only 
to a finite order. In Sec. Ill we investigate the inclusive decay rates of D- and B-mesons 
using the MS scheme. We determine the quark masses and \V c b\ in a manner that is very 
similar to what was done in Refs. |S|-|7| and get a value of \V c b\ that is by about 10% smaller 
than the results of Refs. 0-0. This difference results solely from the transition of the on- 
shell scheme to the MS scheme. Special attention is paid to the correct determination of the 
renormalization scale-dependence. While it is huge in the inclusive rates itself, it mainly 
drops out in \V c b\. Finally, in Sec. IV, we discuss the results and draw some conclusions. 
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II. THE MASS PARAMETER IN HQET 



HQET allows a systematic expansion of QCD observables in inverse powers of a heavy 
quark mass. It is founded on the fact that the long-distance strong interaction of a heavy 
quark with mass mq is independent of its spin and flavour in the limit mq — ► oo. This 
heavy quark symmetry is violated by the finiteness of the heavy quark mass and by short- 
distance QCD interactions. The first correction is taken into account by including terms 
proportional to powers of 1/m.Q, while the effect of hard gluons can be incorporated in a 
perturbation series in a s {mq). Usually these steps are treated independently: At first the 
heavy antiparticle field is integrated out in the path integral and the resulting expression 
is expanded in 1/m.Q yielding a series of local terms in the Lagrangian [IlJ. Then Green 



functions in QCD and HQET are matched at a scale \x ~ mq at which both perturbative 
QCD and HQET are valid. In this step short-distance corrections are taken into account by 
calculating loop corrections to the matching coefficients. Here one has to worry about the 
definition of the mass uiq. The semileptonic inclusive decay rate discussed in Sec. Ill contains 
mq in the fifth power and the proper definition of mq is of significant phenomenological 
importance. 

A mass parameter in a physical observable calculated with the help of HQET can origi- 
nate from two very different sources: Consider the matching of a QCD Green function Gqcd 
to its HQET counterpart Ghqet at some scale fi ~ ttiq: 



GqcdOq,^) = C(m Q ,ii) 



^hqetCaO H ^hqetCaO + O [ — 2" 

m Q V m 6 



(2-13 



In Eq. ( |2.1| ) mq can enter the Wilson coefficient C, which comprises the short-distance 
interactions originating from scales larger than the matching scale \x. The inverse masses 
multiplying Ghq ET , however, is the expansion parameter of HQET multiplying the different 
orders of interaction operators in the HQET Lagrangian. 

The mass rriQ in the Wilson coefficient C, which we will discuss first, is a short-distance 
quantiy and therefore clearly equals the running current mass 7tt.q(//) in the renormalization 
scheme chosen for the calculation of Gqcd- This fact becomes very transparent in the 
example of a QCD Green function whose Feynman graphs contain external heavy quark 
lines and internal lines corresponding to particles with masses mi ^> m 2 ^> . . . mq: We first 
set \i — fii ~ mi and calculate the diagrams in the full standard model and in an effective 
field theory in which the heaviest internal particle is integrated out to obtain some Wilson 
coefficient C(mi,/xi). Whenever we cross a particle threshold in the renormalization group 
evolution of C to lower energies we have to repeat this factorization, which successively puts 
the masses of all internal heavy particles into the Wilson coefficient C. They enter C as 
running masses m^/ij) evaluated at scales //j « m^. When we approach \i f& mq we match 
the Green function in an effective QCD theory to a HQET Green function as displayed in 
( |2.1D , which adds mq to the set of heavy masses appearing as arguments of C. There is no 
point in treating this mass differently from the others in C: The short-distance coefficient 
should not contain any information on the long-distance interaction, i.e. whether the latter 
is encoded in a Green function calculated in HQET or in some other effective theory of 
QCD. 
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In a physical observable calculated with the help of HQET positive powers of mq obvi- 
ously stem from the Wilson coefficient. In spite of this the inclusive decay rate discussed 
in Sec. Ill is usually entirely expressed in terms of the pole mass. The leading term in the 
1/mg-expansion is known to coincide with the parton model result, which has been obtained 
in Jl2| from an older QED calculation. In the latter an on-shell renormalization was used 
for the external fermion so that the result was expressed in terms of the pole mass. In QCD 
this scheme is not very adequate, because quarks are confined and their mass is not directly 
related to any observable. Nevertheless, as long as one stays within the parton model, which 
does not distinguish between the heavy meson and the heavy quark, the on-shell scheme is 
as good as any other scheme, to which one can easily pass by expressing the pole mass in the 
result by the corresponding current mass. By using HQET, however, we want to estimate 
corrections to the parton model result stemming from the binding of the heavy quark in the 
meson, and the correct tool to use is Wilson's operator product expansion (|2.1| ). 



In the case of the inclusive decay rate the standard model diagrams to be calculated for 



the left hand side of ( |2.1|) are the heavy quark self-energy graphs with the light quark and 
the lepton pair in the intermediate state. Its imaginary part determines the desired rate in 
the parton model. In general one would like to include renormalization group improvement 
between the scales fi ~ m^ and fi ~ mq and would match the standard model diagram 
first at fi ~ mw to a corresponding diagram in an effective theory in which the W-boson is 
integrated out, which is then matched at fi ps mq to the self-energy graph in HQET. In the 
case of the semileptonic decay rate, (|3.4|) , however, no such improvement is necessary due 
to the vanishing anomalous dimensions of the corresponding operators, and the standard 
model diagrams are directly matched at /i ~ mq to the HQET Green function (Ai\hTh\A4) 
and to matrix elements of the HQET operators subleading in 1/mq. Here Ai denotes the 
heavy meson, h is the heavy quark field in the effective theory and T is the appropriate Dirac 
structure. In this step one gets the running mass of the heavy quark raised to the fifth power 
evaluated at the scale /i. We stress that one is not forced to choose fi = mq exactly (see 
e.g. Neubert in The choice of fi determines the separation of short and long-distance 
physics, the interaction from scales larger than /i is contained in the Wilson coefficient, which 
contains the running quark mass. In an calculation to all orders the result does not change 
with the variation of /i. In practice one works with a truncated series and the dependence 
of the theoretical prediction on the scale \i estimates the reliability of the calculation. We 
will use this tool extensively in the phenomenolocical analysis of Sec. III. The use of the 
pole mass in the Wilson coefficient corresponds to the matching at an unnaturally low scale 
/i ps 0.5 mq. In the case of the c quark this scale is too small for perturbative QCD to be 
trusted. 

We will now discuss the nature of the mass parameter, whose inverse powers appears in 
the square brackets in (|2.1|) and stems from the HQET Lagrangian. 

The starting point in the derivation of HQET is the decomposition of the heavy quark's 
momentum Pq according to 

P£ = m Q v" + k», (2.2) 

where mq is the heavy quark mass, is the four-velocity of the hadron and k^ 1 is the 
residual momentum. The latter is usually constrained by the condition 
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F = 0(A QC d) - const., (2.3) 

TOq— >oo 

so that finally physical observables X are expanded in a power series in Aqcd/^q : 

X~Ea*(^Y\ (2.4) 
T V m Q / 

In Eq. (|2.3| ) we have emphasized that in HQET the constraint = O(Aqcd) is un- 
derstood as a scaling relation: The residual momentum does not diverge when the limit 
rriQ — > oo is performed. Obviously ( [2.3|) cannot hold for every definition of the mass, be- 
cause a change in the QCD renormalization prescription for rriQ may redefine fc M in ( |2.2| ) by 
an amount involving a s (rriQ)mQV fl , which scales like ttiq/ lnrriQ rather than staying constant. 

Authors addressing the definition of the mass indeed only allow on-shell-like renormal- 
ization conditions for vtlq ||,|,13]. Let us recall their key arguments: The QCD quark 
self-energy 

E = rriQ^i + - m Q )E 2 (2.5) 
can be sandwiched between two projectors = (1 + and be expanded as 
P+EP+ = P+ [mgE, +vkZ 2 ] P+ 

OL s Cf 



An 



m Q E 10 + v-k (En + E 20 ) + O 

\ m Q , 



P+ + O (m Q a 2 s ) . (2.6) 



Here the one-loop contribution to Ej is expanded with respect to v-k/rriQ and a s Cp/(An) Ejj 
is the j-th coefficient. Now E has to match the HQET self-energy v ■ k E up to 0(v ■ k/mo). 
For this an on-shell-like renormalization condition 

Ei| m 2i_p2 =p = (2.7) 

with p pa fixed as mq — > oo was required in ||. Otherwise the term mgE 10 in ( ^.6| ) would 
diverge as rug oo and thereby could not be matched to any HQET Green function, which 
scales at most as a constant due to ( |2.3| ). In the following we set p = 0, so that rriQ in 



equals the pole mass, when the constraint (|2.3| ) is imposed on k^ 



From a phenomenologist's point of view the appearance of a quark pole mass is unsatis- 
factory, because quarks do not exist as free particles, so that the pole mass is no observable. 
Any QCD calculation relates observables to running current masses, which are contained 
in Wilson coefficients. The pole mass has to be calculated in terms of the current mass 
extracted from some experiment. Eq. (|2.2j) implies that the pole mass is calculated to all 
orders in perturbation theory: By truncating the perturbative series defining m pole at order 
n one picks up an error of the order a" +1 (mg)mQ = 0(itiq/ ln n+1 mg) in the residual mo- 
mentum, so that it would not stay constant in the limit ttiq — > oo. Moreover, the authors 
of |10| have proven that the pole mass suffers from an extra infrared renormalon situated 
at u = 1/2 in the Borel plane, so that the result depends on the summation prescription 
chosen for the divergent perturbative series and is ill-defined by a term of order Aqcd- 
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We therefore want to avoid the use of the pole mass in the first step of the construction 
of HQET and use an arbitrary current mass instead. The final result, however, will be the 
same as with the pole mass as the starting point, but an alternative derivation might be 
illustrative. 

If rriQ denotes the current mass in an arbitrary renormalization scheme of QCD, (|2.3|) 
can no more be imposed on /c M and we have instead 



k" ~ 0(a s (m Q )m Q ) = O 



(ln(m Q /A Q cD) 



m Q 



(2.8) 



because a change of the renormalization prescription for uiq redefines k^ by a quantity 
involving a s (rriQ)mQ = 0(rriQ/ lnmg). Nevertheless the heavy quark expansion of QCD 
Green functions remains reasonable, because one still has k^/rriQ — > for rriq — > oo, so that 
v ■ k/rriQ is a small parameter. The heavy quark field h v (x) = P^Q(x) exp(irriQV ■ x) is now 
unambiguously defined, because it involves the current mass. 

We can now investigate how the matching of QCD and HQET has to be performed in 
our approach: Consider the (non-truncated) two-point function G2 in QCD. The tree-level 
part reads 



PtG^P 



iP 



v 9 
p2 



m Q P + 

m 2 n » 



iP+ iP+ 



v ■ k 2m, 



1 - 



(v ■ kf 



O 



1 



m Q In mQ 



In the 0{. . .)-part we have taken into account that fc M scales according to 
part yields for \x = tjiq: 



(2.9) 
D. The 0(a s )- 



Jr., Oro r., 



-iP2 



j + ffl 9j + m Q p 

2 v 



+ 



m 2 Q p 2 



m 



Q 



-%P- 



a s (m Q )C t 
Ait 



(v ■ ky 



J 10 



v ■ k 



j 11 



E 20 + 62 + 1 



(v ■ ky 



no 



+ 



m Q hrniQ 



(2.10) 



where — <5 m Ci?a s /(47r) and S^Cpots/ (Air) are the mass and wave function counterterms in the 
QCD Lagrangian. Eqs. ( |2.9p and ( |2.1(J| ) illustrate that with ( |2.8|) the 1/mQ-expansion now 
yields a series in which each term is still suppressed compared to the preceding one, but 
only by a factor of order 1/ h^mg/AQCD) rather than of order Aqcd/^q- 

The new feature in the matching of (|2.9| ) and ( |2.10| ) to the corresponding expressions in 
HQET is the appearance of the term E 10 — S m in ( |2.1U| ) which equals zero in the on-shell 
scheme. To accomodate for it we need an additional term 



Amh v h v 



in the HQET Lagrangian with 
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Am = ^frriQ (8m - S 10 ) + O (m Q a 2 s ) . (2.11) 

The HQET two-point Green function with one insertion of —Amh v h v then matches the 
first term in fl2.10| ). When working to n-th order in a s , (|2.11|) must be adjusted to cancel 
the mQEi-terms to order a™. Clearly Am in the above one-loop example, Eq. ( |2.1C| ) , is 
nothing but mjjj le — mg. In the MS scheme 

Am(/x) = <ti mQ (\+ In 4 V (2.12) 

Am is a purely short-distance quantity and does not require any definition of m po i e beyond 
perturbation theory. Consequently we treat the "residual mass term" —Amh v h v as an inter- 
action vertex and not as a part of the kinetic Lagrangian. The other term involving Eio — 8 m 
in fl2.10|) stems from the 1/mQ-part of the propagators between which mgSi is sandwiched. 
They clearly match the HQET two-point function with one insertion of —Amh v h v and 
one insertion of the sum of the usual 1/mg-subleading operators 1/ {2mQ)h v {iD) 2 h v and 
—l/(2mQ)h v (iv ■ D) 2 h v . It is evident how the matching of mgS 10 — 5 m to —Am works to 
higher orders in v ■ k/mq. The matching of the wave function counterterms is not different 
from the conventional approach with Am = 0. 

The difference between the residual mass term introduced in [HJ and Am in ( |2.11| ) is 
that the former was constrained to be of order Aqcd while Am is of order a s m,Q. The 
appearance of a positive power of itlq in the HQET Lagrangian has a dramatic consequence 
for the 1/mQ-expansion of Green functions: With the use of the pole mass inverse powers 
of rriQ and logarithms of itiq are neatly separated, so that the 1/mQ-expansion coincides 
term-by-term with an operator product expansion in Aqcd/^q as in Q2.4|) . By using a 
current mass as the expansion parameter and introducing Amg in the HQET Lagrangian 
we have reshuffled the two expansions in I/tuq and « s (mg). In order to arrive directly 
at an operator product expansion (|2.4|) one necessarily has to separate the heavy quark 



mass and self-interaction from the interaction with the light degrees of freedom, which is 



of order Aqcd- The authors of |I(J have shown that this separation is ambiguous beyond 
perturbation theory and can introduce a residual mass term in the HQET Lagrangian which 
is of order Aqcd and is a remnant of the long-distance self-interaction of the heavy quark. 
As for the short-distance contributions to the self-energy of the heavy quark, this separation 
is, however, possible and mandatory to achieve the abovementioned separation of powers 
and logarithms of itlq directly^] 

By using an arbitrary current mass as the expansion parameter we have put the quark 
self-energy in the chosen QCD scheme entirely into the HQET Lagrangian. 

Nevertheless one can easily recover the correct operator product expansion: To this end 
we investigate how the hadronic parameters of HQET transform under a change of the 
renormalization scheme used to define tuq. 

Consider two schemes differing by a finite renormalization of the mass: 



We thank M. Beneke for clarifying this point. 
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m' Q = z m m Q = m Q + Sm Q 



Am' Q = ArriQ - 5m Q . (2.13) 



The choice Amg = Sttiq transforms to the pole mass. The transformation ( |2.13| ) modifies 
the HQET quark field as 

h' v (x) = h v e ivxSm Q. (2.14) 

The parameter A measures the mass difference between meson and quark. One easily finds: 

- ivd(OffTh' v \M(v)) 
(0\q'Th' v \M(v)) ' 

= A-5m Q , (2.15) 

so that the meson mass 

M Q = rriq + A + ( ^] (2.16) 
V m Q J 

is scheme-independent. In an arbitrary scheme A contains the quark self-energy, whose 
short-distance part depends on tuq, so that we have to add a flavor label to A and write in 
the e.g. MS scheme for \x = tuq: 

Af = Apoie + \^^m Q + O , (2.17) 

where A po i e is the usual A-parameter which corresponds to the choice of ttiq = m po \ e and 
obeys the heavy quark symmetry. So far we have just trivially pushed self-energy contri- 
butions from the mass into Aq. Reference to the pole mass is for the first time made in 
( p.!7| ), which expresses that only with an on-shell renormalization for tuq the quantity A 
is of order Aqcd- From |10[ we know that the separation of m po i e and A po i e is ambiguous 
by long-distance terms of order Aqcd- With our formalism we clearly do not resolve this 
ambiguity, but just make clear that it does not affect a practical phenomenological analysis, 
where one would like to use ( |2.16| ) to extract, say, A pole from an analysis of the D-meson 
system and to insert it into a prediction for the B-meson system. In the first step one ex- 
tracts the current charm quark mass from some observable which is calculated with radiative 
corrections to some order a™. Clearly in Q2.17Q one will only subtract the quark self-energy 
to the same order, so that the prediction for the calculated observable will then have an 
error of the order One is never faced with the problem to define the pole mass beyond 

perturbation theory. Even if we did the calculation to a very high order, so that we saw the 
renormalon-induced divergence of the perturbative series for the self-energy in Eq. ( p.!7|) , 



we would naively expect that the renormalon ambiguity cancels in the prediction, if we fixed 
the summation prescription for the perturbation series in both the D- and B- analysis in the 
same way. 

We will next discuss the scheme-dependence of the parameters Ai and A2, which appear 
at order I /nig and parametrize the matrix elements with two heavy meson states of equal 



S 



velocity. At this order the residual mass Am has to be taken into account. We first want 
to define Ai and A2 in a scheme-invariant way, i.e. such that they are independent of Am. 
From WM we know that the latter can be obtained by replacing the covariant derivative 



by .D M + iv^Am. Indeed, for a heavy meson Ai with mass Mq we find 

2M Q \ 1 = (M\h v (iD -vAm) 2 h v \M)(l + (D(l/m Q )), (2.18) 

6M Q A 2 (//) = ( M I h v 9 - a, u F^h v \M)(1 + 0(l/m Q )) 



to be invariant with respect to the transformations ( p. 13 ) and ( 2.14]) . Consider now some 
observable whose 1/mg-part looks for itlq = m po \ e like 

Ax 



^■polc 



(a + ba 8 + O (cty) (2.19) 



with some constants a and b plus a similar term with X 2 . In an arbitrary scheme one 
gets the current mass instead in the denominator of Q2.19p , but the second order in the 



1/mg-expansion involves matrix elements with one insertion of Am and two insertions of 
the 1/mg-operators. Since Am/mg = 0(a s ), these contributions are of the same order as 
the first order radiative corrections in (|2.19|) . Renormalization scheme-invariance of QCD 
requires that both terms combine to a scheme-invariant quantity: 

1 A "' ■ " 0(ai)). (2.20) 



m Q m 2 Q m W e 

If we want to include explicit radiative corrections to order n in (|2.19|) , we must calculate 
Am to order a™ and have to take into account matrix elements with multiple insertions of 
Am up to order l/mg +1 , so that the reordering of the expansions in 1/mg and a s takes 
place over n orders in the 1/mg-expansion. Hence the final result is the same as with the 
use of the pole mass from the very beginning. Yet in our derivation the pole mass appears 
explicitly only in the last step of the derivation. One has to calculate it only to the same 
(finite) order in perturbation theory as the explicit radiative corrections to the corresponding 
term in the observable (|2.19| ) under consideration. 

In the sketch given above the ambiguously defined pole mass to all orders does not appear 
in the definition of the HQET field h v , instead the self-energy —Am of the quark is encoded 
order by order in a s in the HQET Lagrangian. We have made plausible that the use of 
HQET in a phenomenological analysis, which includes radiative corrections to some finite 
order to the various terms in the 1/mg-expansion, does not require the calculation of m po i e 
beyond the same order in its perturbative series. 

Let us emphasize again that our formalism of HQET with the current mass as the 
starting point leads exactly to the same result as the conventional one, if one uses the 
pole mass calculated perturbatively to the same order as the explicit radiative corrections 
multiplying the 1 / mg-term under consideration. Our statements concerning the role of the 
mass appearing in the Wilson coefficients as, e.g., the mg-term in the inclusive decay rate, 
lead, however, to numerical effects in phenomenological analyses as demonstrated in Sec. III. 



9 



III. SEMILEPTONIC INCLUSIVE DECAYS 



Let us now turn to an important application of the considerations of the last section. 
Quite recently, the heavy quark expansion has been extensively applied to the problem of 
inclusive heavy meson decays, both semi- and nonleptonic ones ||15|| . Two major results 



were obtained, namely first that inclusive decays are determined by the free quark decay in 
leading order in the heavy quark expansion, and second that non-perturbative corrections 
to the free quark decay picture are suppressed by terms of order l/m^; there are no terms 
of order 1/itlq. These results allow an immediate application to the determination of the 
CKM matrix element \V cb \. This task was tackled in the papers Ref. where the following 
procedure was adopted: first the c quark mass m c is determined from the branching fraction 
B(D — ► Xev) in dependence on the unknown non-perturbative parameters. Then HQET 
is invoked in order to fix the b quark mass m b from the known value of m c . Insertion in 
the theoretical expression for B(B — > X c eu) finally yields \V^,\ [||||, or, assuming \V c b\ as 
known from other sources, bounds on the non-perturbative parameters and quark masses, 
respectively, cf. [[/]]. In the present analysis we will closely follow the sketched method, but 
work in the MS scheme and keep track of possible scheme-dependences. 

The non-perturbative corrections to the free quark decay picture can be expressed in 
terms of two matrix elements, 

2m B X 1 = ( B | b v {iD - vAm)\ \B)(1 + 0(l/m b )), 
6m B X 2 (fi b ) = ( B | 6«| a^F^b* \B)(1 + 0(l/m b )), (3.1) 

where b v is defined as b v = e imbVX b(x), b(x) being the b quark field in full QCD, and is the 
four-velocity of the B meson. The corrections of order 1 / m b on the right-hand side account 
for the fact that physical states in QCD and HQET differ by terms of order I/tuq. The 
above quantities have an immediate physical interpretation by virtue of the scheme-invariant 
relation 



Ai + 3A 2 (/x b ) ^ / as 

2m b \m b m b 



m B = mb + A b - ^\- 2 ^ +0[^,- 2 \ (3.2) 



where \i/(2m b ) plays the role of the kinetic energy of the heavy quark's Fermi motion 
inside the meson and the term in X 2 , which has non-zero anomalous dimension in HQET, 
accounts for its spin-energy in the chromomagnetic field. fi b is the matching scale of HQET 
onto QCD. A b can be interpreted as the binding energy of the light degrees of freedom in 
the meson. From ( |3.2|) one readily infers 

A 2 (/i 6 ) w -{m 2 B * - m 2 B ) w 0.12 GeV 2 (3.3) 

which is true in any renormalization scheme. Unfortunately, A i is no observable, but has to 
be determined within some model-calculation. Recently, the value Ai = —(0.6 ± 0.1) GeV 2 



was obtained from QCD sum rules [16]. In view of the criticisms raised in |17| , we neverthe- 
less prefer to leave Ai as an open parameter to be varied within the interval [—0.7, 0] GeV 2 
where the upper bound is taken from ||, the lower one from the sum rule determination 
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The heavy quark expansion of the decay rate of the semileptonic decay B — > Xeu with 
unobserved X in the final stated then reads : 



1927T 3 



1 2a;' J |ffl c 



3 7T 



9A 2 / m c 
/2 —7? 



2(m^) 2 



+ 



O 



Ai 



2(mF) 2 



/i 



n? 



mi 



la 



R\2 



mi 



m 



(3.4) 



Here we have marked all quantities depending on the renormalization prescription by the 
superscript "R" . j\ and f 2 are phase-space factors given by 



AO*) 
h{x) 



lx 2 + 8x 6 - x 8 - 24x 4 lnx, 



8 5 
1 - - x 2 + 8x 4 - 8x 6 + - x 8 + 8x 4 \nx. 
3 3 



(3.5) 



_,on— shell 



was first analytically calculated in Ref. [IS 



h{x)g 



on-shell , 



X 



25 239 



25 



4 



+ 4x 4 (36 + x 4 )lnx + (1 - x 4 



1 - an — x A + — x 4 + 2x^ 20 + 90x^ --i 4 + -/ lni 
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17 64 



x 



— x 4 ^ln(l — x z 



3 3 

- 8(1 + 30x 4 + x 8 ) In x ln(l - x 2 ) - (1 + 16x 4 + x 8 ){6Li 2 (x 2 ) - vr 2 } 

1 



32x 3 (l + x 2 Kvr 2 -4Li 2 (x) + 4Li 2 (-x) -4 In a; In 



x 



1 + x 



(3.6) 



g Mb can be obtained from ( |3.4j ) by expressing all scheme-dependent quantities in the MS 



scheme. To first order in a s , we thus only need the relation 

,2 \ 



m P° lc = 7716 (fl) |l 

Insertion in (|3.4j ) yields 



O's 



4 /i 
+ In — 



OK 



(3.7) 



S MS (*) 



<T^(x) + 3x \nx dlnf / X) - 10 + ^ In 4- 

ax 2 /r 



(3.8) 



Numerically, we find the values given in Tab. |. From the comparison of g R in both schemes, 
it becomes evident that higher order corrections in a s are of paramount importance if one 
is willing to obtain reliable predictions. The effect of higher order terms is conventionally 



2 In the following we identify B — > Xev with B — > X c eu where the b quark decays into a c quark 
since charmless decays B — > X u ev can safely be neglected because of \V u b/V c t,\ 2 < 1% [18]. 
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estimated by allowing the renormalization scale /i to vary within, say, ra&/2 < // < 2m},. 
Let us compare this with the possible improvement achievable by a future calculation of the 
term in a 2 s . It is of the generic form 



a 



^| Cl ln 2 ^ + c 2 ln^ + C3 



mi 



(3.9) 



where both ci and c 2 are completely determined by the terms of lower order in a s and only 
C3 remains to be calculated. Expanding the anomalous dimension of the running MS mass 
as 



7 



and the QCD /3-function as 



(3.10) 



(3.11) 



we find^ 



Cl 



<~2 



128 



7o m (5 7o m -2/3 ), 



fj,=m b 



(3.12) 



In Fig. |l] we illustrate the resulting scale-dependence of the branching ratio T(D — > 
Xev) for 1 GeV < fi < 2 GeV. We show both the next-to-leading order (NLO) result 
( j3.4|) and the NNLO result with the a 2 s terms (|3.9|) included where the unknown constant 
C3 is varied in the interval [—10,10]. As can be read off the figure, the scale-variation 
in the NLO result reproduces approximately the inherent uncertainty from the unknown 
constant C3. Still, the remaining scale-uncertainty is noticeable and limits the accuracy 
achievable in the determination of m c (m c ) from T(D — > Xev). We illustrate that point 
in Fig. ^ where m c (m c ) as determined from T(D —>■ Xev) via Eq. (3~3) is plotted as a 
function of the renormalization scale \i and for different values of the input parameters 
A: (Fig. 1(a)), B(D - Xev) (Fig. 1(b)), m s (lGeV) (Fig. 1(c)), and Ag, (Fig. 1(d)). 
In Fig. ^(e) we also show the possible effect of 1/mg corrections estimated by adding to 
and subtracting from (|3.4j) the 1/ttiq terms to the power 3/2. As input parameters we 
use0 B(D -> Xev) = 0.172 ± 0.019, r D+ = 1.066 ps, \V CS \ = 0.9743, \V cd \ = 0.221 0, 
GeV) = (0.20 ± 0.05) GeV, cf. e.g. @], and = (300 ± 50) MeV @. In contrast 



3 The numerical values of the coefficients 7™ and /3, can be found in any good textbook on QCD. 

4 The experimental errors on t d +, \V cs \, and \V c d\ are so small, that the combined error of all 
experimental quantities is completely determined by the error of the branching ratio, and thus can 
safely be neglected. 
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to the analysis done in 0] , we have varied \i in the range 1 GeV < < 2 GeV only, since 
we believe that the perturbative expansion becomes highly unreliable at smaller scales. We 
cannot follow the arguments of Luke and Savage that /i could be associated with a "typical 
energy" of the emitted lepton pair and be as small as m c /3 ~ 0.5 GeV! In our approach, 
H has to be identified with the matching scale of HQET to QCD and cannot be attributed 
any physical meaning. The only (loose) condition to be imposed on [i is that it should be of 
order m c . We would also like to mention that we are rather suspicious about how the scale- 
uncertainty is handled in Refs. [0-0] since in these papers only the scale in a s (fi) is varied 
(and the MS scheme used in doing that!), but the quark mass kept fixed. This corresponds 
to a scheme where the running of the quark mass is cut (7™ is put zero in ( |3.12| )), which is 
a permissible scheme, but does not seem a very reasonable one. At least it is not the correct 
way how the QCD scale-uncertainty should be estimated in the on-shell scheme. 

From Fig. |2| we find that m c (m c ) is clearly most sensitive to the choice of /1. In particular, 
we note that a variation of Ai within the conservative range [-0.7, 0] GeV 2 affects the value 
of m c (m c ) only as much as a variation of B(D — > Xev) within one standard deviation 
and nearly as much as a variation of the value of the strange quark mass by 50 MeV. 
These observations seem to cast some doubts on the procedure employed in where the 
determination of m c was used to derive bounds on Ai and A. Taking all together, the heavy 
quark expansion yields 

m jp(m c ) = (1.35 ± 0.20) GeV for 1 GeV < /1 < 2 GeV. (3.13) 

We next use the information gained from the decay of the D meson to determine m b . To 
that end, we use Eq. (|3.2|) and the relation 



- 4 a s (m c ) j- 4a s (m b ) 

A c -~ m c (m c )=A b -- m b (m b ) (3.14) 

o 7T on 

which is valid in the MS scheme and can be obtained from Eq. ( [2. 17j ). The results are 
displayed in Fig. 0, both for m b as function of Ai with fixed scale \i = m c , Fig. |3](a), and for 
fixed Ai = GeV 2 as function of the scale /x/m c , Fig. 0(b). Again we demand /i > 1 GeV, 
which yields /i/m c > 0.8. The dashed lines give the uncertainty of m b due to the experimental 
error in B(D —* Xev). Once more we observe a strong dependence on the renormalization 
scale /1, which is however milder than for m c due to the higher scales involved. In contrast 
to the c quark mass there is also a sizable dependence of m b on Ai which is of the same order 
as the /^-dependence and originates in the relation ( j3T2[ ). We obtain 



m b (m b ) = (4.6 ± 0.1 ± 0.1 ± 0.1) GeV, (3.15) 

the first error being due to the dependence on Ai, the second one due to /i-dependence, and 
the third one combines the dependence on the s quark mass, the branching ratio B(D — > 
XHandAg, 

In determining V cb , we in addition need the experimental branching ratio for B — > X c ev, 
B(B — > X c ev) = (10.7 ± 0.5)%, quoted in |fI8|| , and the B lifetime tb, where we use the 



most recent value t b = 1.49 ps quoted in [OT. In view of the changes the value of tb has 
experienced in recent years, cf. ||18|| , we prefer not to include its experimental uncertainty in 



the error analysis, but to give our results for the above fixed value of t b . In Fig. we show 
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B(B — > X c ev) as function of /x with = 0.04, m c {m c ) = 1.35 GeV, mbirrib) = 4.5 GeV, 
Ai = GeV 2 and A 2 (/i&) = 0.12 GeV 2 . Analogously to Fig. [I], the solid line obtained from 
Eq. ( p.4|) includes 0(a s ) terms with varying /x, whereas the dashed lines also include the 
0(a 2 s ) corrections ( p.9|) with c 3 varied within [—10, 10]. Compared with Fig. [I], there is still 
a scale-dependence visible, but it is weaker due to the higher scales involved. Still, the 
possible effect of a future calculation of C3 can reliably be mimiced by varying /x within 
m b /2 < /x < 2m b . 

We are now in a position to determine |V^,|. Inserting correlated values of m c and mj, 
into Eq. ( |3.4j ) and comparing with the experimental branching ratio B(B — > X c ez/) we 
obtain the values shown in Fig. |5]. An analysis of the different sources of uncertainties 
shows that they are dominated by the dependence on Ai, which proves the strongest one, 
and on B(B — > X c ev). In Fig. |(a) we thus plot \V c b\ as function of Ai for fixed scales 
fib/ m b = Hc/ m c = 1 and in Fig. |5|(b) as function of the scale (jl/itiq for fixed Ai = OGeV 2 ; 
we also show the values of \V c b\ resulting from a variation of B{B — > X c ev) within one 
standard deviation (dashed lines). For the convenience of the reader, the numerical values 
are given explicitly in Tab. [TT]. The strong dependence of \V c b\ on Ai is entirely due to the 
behavior of m& as function of Ai since the explicit Ai/m 2 terms present in Eq. (|3.4j ) give only 
tiny contributions. The rather marginal scale-dependence of \V c b\ visible in Fig. |5|(b) can be 
explained in the following way. In Tab. |TTT| we give for different scales /x/ m Q the values of 
the quark masses, their ratio that enters the phase-space factor, the value of the latter one, 
the branching ratio divided by phase-space and the branching ratio itself, calculated with 
\V c b\ = 0.04. Although the quark masses differ rather drastically, the ratio m c (/ift) / 'm&(/x&) 
is not very sensitive to /x^ and also the phase-space factor fi varies by only about 10% in 
the range 0.8 < fib/m b < 1.8. The branching ratio divided by phase-space as well as the 
branching ratio itself are even less sensitive to the scale. We thus find that the change in 
rribinb) is nearly completely compensated by a corresponding change in the terms in a s (/Xfe) 
which is far from being trivial. From the figures we read off 

/ \ 1/2 

\V cb \ TB = 0.036 ±0.002 ±0.001 ±0.002, (3.16) 

\1.49 ps J 

where the first error includes the dependence on Ai in the interval [—0.7, 0] GeV 2 , the second 
one the scale-uncertainty, and the third one all other uncertainties in the parameters, i.e. 
in the branching ratios, m Sl and In principle, there is also an uncertainty due to the 

scheme-dependence of our analysis which is however difficult to estimate. If the on-shell 
scheme was trustworthy, we could take the difference of our results and those of as 
rough estimate for that "hidden" theoretical error. For lack of any other calculation in a 
reasonable scheme we are however forced to leave that point open to future investigations. 
The above value of jV^I has to be compared with 

/ \ 1/2 

lVcbl [m^) = (0 - 046 ± °- 008) [5] ' 
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tb 



1/2 



c6 



1.49 ps 



0.037 ±0.007 



(3.17) 



The value quoted in [21 is the most recent one obtained from the spectrum of the exclusive 
decay B — > D*ev fitted to the Isgur-Wise function using different shapes. We observe that 
our result nearly coincides with the one from exclusive decays, but is smaller than the ones 
obtained in where the on-shell scheme was used. Let us close this section with a short 
comment on the prospects of a future reduction of the errors in ( 3.16Q . To begin with, 
experimentalists are challenged to carry out more accurate measurements of the branching 
ratios, especially of B(B — > X c eu). Combining with more sophisticated determinations of 
22| and A^ will reduce the third error. As for the second error due to scale-dependence 



m* 



there is no simple remedy at hand and we hardly can imagine any reliable method to be 
invented in the near future to fix the scale except for a calculation of the 0(a 2 s ) corrections 



to Eq. fl3.4| ) which is a truly formidable task. So the main efforts should be concentrated 
on the determination of Ai. Although in there was proposed a method of relating Ai to 
moments of the so-called "shape-function" determining the end-point region of the lepton- 
and photon-spectrum in the decays B — > X u ev and B — > X s j, respectively, at present 
that approach suffers from missing both measurements and any formulas going beyond 
0{oPg). Yet the introduction of the shape-function leads to the constraint A < GeV 2 which 
is in agreement with the results obtained in [16] and [23|, but in contrast to the earlier 
24"f , where a positive value was obtained. The status of the careful analysis 
with QCD sum rules is unfortunately not beyond any doubt, cf. fl7|| , so a 
clarification of that point would be welcome. If the result of [|ll]] could be confirmed, fl3.16p 
would read [^1(^/1.49 ps) 1 / 2 = 0.034 ± 0.001 ± 0.002. 



determination 
done in 
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IV. SUMMARY AND CONCLUSIONS 

In applications of the heavy quark expansion of QCD to physical quantities involving 
positive powers of the heavy quark mass the proper definition of the mass parameter is of 
supreme phenomenological importance. We have made clear that the heavy quark mass 
parameter stemming from the Wilson coefficient in the matching of perturbative QCD and 
HQET Green functions, in a rigorous QCD context enters the calculation as the running 
current mass rather than the pole mass. It is evaluated at the scale at which we perform 
the matching. This scale must be of order of the heavy quark mass, but can be varied to 
estimate the error caused by the truncation of the perturbative series. The choice of the 
on-shell scheme, which is widely used also for the mass parameters in the Wilson coefficient, 
corresponds to a special choice for /i, which is unnaturally low in the case of the c quark. It 
is well-known that the heavy quark expansion reproduces in leading order the parton model. 
With the proper identification of the mass parameter stemming from the Wilson coefficient 
we also recover the familiar renormalization scheme- and scale-dependence of the parton 
model. 

The natural expansion parameter of the HQET Green functions, however, is the pole 
mass, which was recently found to be ambiguous when defined beyond finite orders in pertur- 
bation theory. Since HQET attempts to include non-perturbative effects, a phenomenologist 
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using HQET has to worry whether this ambiguity may be relevant for his physical predic- 
tions. We have argued that a calculation which includes only radiative corrections to a finite 
order to a given term in the 1/mQ-expansion, requires only the calculation of the perturba- 
tive pole mass to the same order in a s . This fact is plausible, but not completely obvious, 
because the usual derivation of HQET requires the all-orders pole mass in the first step to 
define the heavy quark field h v . 

The proper identification of the mass parameters becomes probably most pronounced 
in the case of semileptonic inclusive decays. Their decay rate depends on the quark mass 
in its fifth power and is most sensitive to any change in the definition of that parameter. 
Here it likewise becomes rather obvious that the on-shell scheme is, even apart from the 
renormalon problem, a rather unnatural one, since it is always the current mass that appears 
in actual calculations, whereas the pole mass has to be put in afterwards by hand, adding 
higher order terms in a s in an uncontrollable manner. We have carefully studied the effect 
of changing the scheme and scale and found that neither the c nor the b quark mass can be 
determined reliably from inclusive decays. In extracting \ V c b\, however, the scale-dependence 
mostly cancels and we found a reduction of the value of \V c b\ by more than 10% compared 
with results obtained in the on-shell scheme By a more consistent error analysis, in 

particular of the scale-dependence, we also could reduce the absolute error of |V^|. We 
also clarified the meaning of the renormalization scale fi entering the parton model and its 
non-perturbative corrections and identified it with the scale at which HQET is matched to 
QCD, a point that was not paid proper attention to in the previous analyses |5|-|7|. 

Finally, we suppose that our considerations could also give some hint at the solution of 



the problem of the "baffling semileptonic branching ratio of B mesons" raised in [5S| and 
that it could be worthwile to attack that problem starting from a different choice of the 
renormalization scheme. 
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FIG. 1. The branching ratio B(D — > Xev) as function of the renormalization scale \i for 
m c {m c ) = 1.35 GeV, Ai = OGeV 2 , \2{rn c ) = 0.1 GeV 2 . Solid line: B = t^T calculated according 
to Eq. (|3.4| ), dashed lines: B according to Eq. fl3.4|) with inclusion of 0(a 2 (/x)) corrections according 
to Eq. ( |3.9[ ). The different dashed lines correspond to three different values of C3: 10, 0, —10 (from 
above). 
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Fig. 2(c) 
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Fig. 2(e) 
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FIG. 2. m c (m c ) as function of the renormalization scale \jl for different values of the in- 
put parameters. The solid lines are calculated with Ai = OGeV 2 , B(D — > Xev) = 0.172, 
m s (l GeV) = 0.2 GeV, = 300 MeV. The (long-, short-) dashed lines are obtained by replacing 
the respective parameter by (a) Ai = (-0.35, -0.7) GeV 2 , (b) B(D -» Xev) = (0.191,0.153), (c) 
m s (lGeV) = (0.15, 0.25) GeV, (d) = (250, 350) MeV. In (e) we have included l/m 3 c correc- 
tions as described in the text. 
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Fig. 4(a) 



Fig. 4(b) 
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FIG. 3. (a): m&(m&) as function of Ai for fj. 
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The solid line is calculated with the same 



set of input parameters as in Fig. The dashed lines correspond to a variation of B(D — > Xev) 
within the experimental error: B = 0.153 (short dashes), B = 0.192 (long dashes), (b): m b {m,h) as 
function of ///m c for Ai = OGeV 2 . The parameter sets used for the curves are the same as in (a). 




FIG. 4. The branching ratio B(B — > X c eu)(0.04:/V c b) 2 as function of the renormalization scale 
H for m c (m c ) = 1.35 GeV, m b (m b ) = 4.5 GeV, Ai = OGeV 2 , \ 2 {m b ) = 0.12 GeV 2 . Solid line: 
B = tbT calculated according to Eq. (EO), dashed lines: B according to Eq. (|3.4D and with 



inclusion of 0{a 2 s (pb)) corrections according to Eq. (3.9). The different dashed lines correspond to 
three different values of C3: 10, 0, —10 (from above). 
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1. The solid line is calculated 



300 MeV. 



(a): I V^fo I (ts/1.49 ps) 1//2 as function of Ai for /j/mg 
> X c ev) = 0.107, B{D -» Xev) = 0.172, m s (lGeV) = 0.2 GeV and 
The dashed lines correspond to a variation of B{B — ► X c eu) within the experimental error: 
B = 0.102 (short dashes), B = 0.112 (long dashes), (b): iV^r^/l^ps) 1 / 2 as function of the 
ratio [i/rriQ for Ai = OGeV 2 . The parameter sets used for the curves are the same as in (a). 
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TABLES 



X 


0.0 


0.1 


0.2 


0.3 


0.4 


0.5 


0.6 


0.7 


0.8 


0.9 


1.0 


g on shell 


3.62 


3.25 


2.84 


2.51 


2.23 


2.01 


1.83 


1.70 


1.59 


1.53 


1.50 


9 MS 


-6.38 


-7.79 


-9.94 


-12.1 


-14.2 


-16.0 


-17.8 


-19.4 


-20.9 


-22.2 


-23.5 



TABLE I. Radiative corrections to the semi-leptonic free quark decay, Eq. ( |3.4|) , in the on-shell 
and in the MS scheme. 



Ai 


[GeV 2 ] 


-0.7 


-0.6 


-0.5 


-0.4 


-0.3 


-0.2 


-0.1 





B 


= 0.102 


0.0329 


0.0334 


0.0339 


0.0344 


0.0350 


0.0356 


0.0362 


0.0367 


B 


= 0.107 


0.0337 


0.0342 


0.0347 


0.0353 


0.0358 


0.0364 


0.0370 


0.0376 


B 


= 0.112 


0.0344 


0.0350 


0.0355 


0.0361 


0.0367 


0.0373 


0.0379 


0.0385 



TABLE II. |T4b|(Ts/L49ps) 1 / 2 for several values of Ai and the branching ratio B(B — > X c ev). 



The same input parameters are used as in Fig. ||(a). 



/V m Q 


m c (m c ) [GeV] 


m b (m b ) [GeV] 


m c (n b )/m b (^L b ) 


h 


Tr B /h 


Tt b 


0.8 


1.28 


4.41 


0.21 


0.72 


0.171 


0.123 


1.3 


1.43 


4.57 


0.23 


0.67 


0.182 


0.122 


1.8 


1.51 


4.65 


0.24 


0.64 


0.179 


0.116 



TABLE III. Quark masses, their ratio, phase space h an d branching ratio Ttb for different 
values of the scale n/mQ. In the last two columns we have put \V cb \ = 0.04. 
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